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ABSTRACT 

We present multi-epoch radio and optical observations of the M7 dwaif 2MASS J13142039+1320011. We 
detect a ~ 1 mJy source at 1.43, 4.86, 8.46 and 22.5 GHz, making it the most luminous radio emission over the 
widest frequency range detected from an ultracool dwarf to date. A 10 hr VLA observation reveals that the radio 
emission varies sinusoidally with a period of 3.89 ±0.05 hr, and an amplitude of « 30% at 4. 86 GHz and « 20% 
at 8.46 GHz. The periodicity is also seen in circular polarization, where at 4.86 GHz the polarization reverses 
helicity from left- to right-handed in phase with the total intensity. An archival detection in the FIRST survey 
indicates that the radio emission has been stable for at least a decade. We also detect periodic photometric 
variability in several optical filters with a period of 3.79 hr, and measure a rotation velocity of vsin/ = 45 ± 5 
km s"', in good agreement with the radio and optical periods. The period and rotation velocity allow us to 
place a lower limit on the radius of the source of > 0.12 R0, about 30% larger than theoretical expectations. 
The properties of the radio emission can be explained with a simple model of a magnetic dipole mis-aligned 
relative to the stellar rotation axis, with the sinusoidal variations and helicity reversal due to the rotation of 
the magnetic poles relative to our line of sight. The long-term stability of the radio emission indicates that the 
magnetic field (and hence the dynamo) is stable on a much longer timescale than the convective turn-over time 
of ^ 0.2 yr. If the radio emission is due to the electron cyclotron maser process, the inferred magnetic field 
strength reaches at least 8 kG. 
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1. introduction 

An unexpected result from radio observations of low mass 
stars and brown dwarfs was the recent disco very of peri- 
odic r adio emission tied to the stellar rotation dBerger et alj 
I200I iHaiiiiian et all llOOel 120071: iBereer et all 12009V Pre- 
viously thought to be incapable of generating and dissipat- 
ing strong and stable magnetic fields, observations in recent 
years indicate that at least 10% of fully convective late-M 
and L dwarfs possess stable and large-s cale magnetic fields 
(iBerged I2OOI 120061: iBerger et alJ l20Toh. W hen combined 
with X-ray and Ha o bservations dMoh a ntv & Basri 2003; 
iReiners & Basril 120081: iBerger et alj |20IO|), Zee man'broad- 
ening of FeH molecular lines in Stokes / (iReiners & Basril 
l2007h . and spe ctropolarimetry in Stokes V (Zeeman Doppler 
Imaging: ZDI: ' Donati et al.ll2006h . a picture of the magnetic 
topology of fully convective ultracool dwarfs is now begin- 
ning to emerge. 

Zeeman broadening observations provide a measure the 
average surface magnetic field, Bf (where B is the mag- 
netic field and / is the covering fraction) and have been 
used to study modera tely rotating objects as cool as M9 
( IReiners & Basril [20101) . These observations uncovered av- 
erage magnetic field strengths of up to a few kG. The ZDI 
technique has been used to study the large scale structure 
and topology of magnetic fields (as it is sensitive only to 
the net polarization in Stokes V) in objects down to spec- 
tral type M8. At the fully convective boundary (spectral 
type of ^ M3) these studies hint at a transition from non- 
axisymmetric toroidal fields to predominantly poloidal, ax- 
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isymm etric field topologies (iDonati et al.ll2008t iMorin et aP 
120081), although the systematic uncertainties and general de- 
tectability of field topologies across a wide range in stellar 
parameters are not fully understood. Similarly, there seems to 
be an increase in the ratio of magnetic fluxes measured in cir- 
cularly polarized to unpolarized spectral lines at this bound- 
ary, apparently confirming that the fully convective objects 
have increasingly large-scale, low-multipole field topologies. 
However, a number of fully convective mid-M dwarfs, poten- 
tially belonging to a younger p opulation, still app ear to pos- 
sess small-scale complex fields (iMorin et al.ll2OI0|). alth ough 
recent observations seem to question this idea (Reiners 2010). 

In the context of magnetic activity, the correlation be- 
tween X-ray and radio emission, which holds across sev- 
eral orders of magni tude in each quantity and over a wide 
spect ral type range dGuedel & Benzl 1 1 9931: iBenz & Guedell 
I1994I), sharply breaks down around spectral type M7-M8 
(Berger etal. 2001; Berger 2002; Berger etal. 2010). X-ray 
and Ha activity also decline precipitously around the same 
spectral type (.Mohantv & BasrilllOOl IReiners & Basrill2008l: 
iBerger et aP l2010l) . Finally, the Ha and X-ray rotation- 
activity relations also break down near the M/L transition 
(Mohanty & Basri 2003; Berger et al. 2008; Reiners & BasrJ 
I2OIOI) . The origin of these changes remains unclear, in par- 
ticular whether they are related to a change in the magnetic 
dynamo itself, or to a change in the way magnetic energy is 
released in the cool and largely neutral atmospheres. 

Radio activity in ultracool dwarfs, on the other hand, ap- 
pears to be equally strong (or stronger when normal ized to 
the bolometric luminosity) th an in early-M dwarfs ( Berged 
l2002ll2006l:lBerger et al.ll2010l) . Among the active population 
of ultracool dwarfs, the handful of objects with periodic ra- 
dio emission exhibit diverse behavior. In particular, the wide 
variation in the temporal behavior, emission bandwidth, and 



fractional polarization suggest that there are a number of dif- 
ferent processes at w ork. Three object s (TVLM 5 13-46546 : 
lHallinanetal] l2006': LS R 1835+3259: iHallinan eFdl 120071: 
and 2M 0746+2000: Ber ger et al.ll2009h have been observed 
to emit short-duration, highly polarized pulses with appar- 
ently narrow bandwidth. On the other hand, 2M 0036+1821 
exhibits periodic variability with a different s tructure, ris- 
ing an d falling in a roughly sinusoidal manner dBerger et alj 
l2005h . In all four objects, the observed periods, which range 
from two to three hours, are in close agreement with the 
known rotation velocities (vsin/). 

Two primary emission mechanisms have been proposed 
to explain the observed radio emission: gyrosynchrotron 
radiati on and the e l ectron cyclotron maser (ECM) insta- 
bility ('Berger ■2006; iTreumannI 120061: iHallinan et all 120081: 
iBerge r et al. 2009). The short-duration periodic pulses have 
been attributed to a coherent process such as ECM due their 
small bandwidth and high degree of circular polarization 
(IHallinan et al.l2 007: Berger et al. 2009). However, the mech- 
anism for the sinusoidal radio emission from 2M 0036+1821, 
as well as the quiescent emission from several ultracool 
dwarfs, is still under debate, partly because of the low levels 
of circular polarization and the much wider bandwidth. 

Here we present multi-wavelength observations of 
a newly discovered radio-emitting M7 dwarf, 2MASS 
J13 142039+1320011 (hereafter, 2M 1314+1320), which 
exhibits unique emission properties among the population 
of radio active ultracool dwarfs. It has the brightest radio 
emission over the widest frequency range of any ultracool 
dwarf to date (detected from 1.4 to 22.5 GHz with a nearly 
flat spectrum), and presents the first clear example of si- 
nusoidal emission (in phase at 4.96 and 8.46 GHz). The 
periodicity is also detected in circular polarization, including 
a helicity reversal that provides insight into the topology 
of the emission regions. Photometric observations reveal 
periodic variability with the same period as in the radio 
band, indicative of a large-scale magnetic spot structure. The 
unique radio properties of 2M 1314+1320, and its location 
at the cusp of the transition in magnetic field properties 
(spectral type M7), make it an invaluable laboratory for fully 
convective dynamos. 



2. OBSERVATIONS 

2.1. Target Properties 

We identified the M7 ultracool dwaif 2M 1314+1320 as 
a radio-emitting source in 2007 during a large Very Large 
Array (VLA^) survey of nearby mid- and late-M dwarfs 
dMcLean & Bergen 1201 lb . It is a high proper motion 
object ((5RA fti -244 mas yr~', (5Dec k, -186 mas yr"'; 
iLepine & Sharal 12005!) . consistent with the you ng disk pop- 
ulatio n, and has a parallactic distance of 16.4 pc dLepine et al.l 
120091) ■ Observations in the / and z bands with the Lucky 
Imaging technique have demonstrated that it is a resolved bi- 
nary with a separation of 0.1" (about 1.6 AU), and a com - 
panion that is fainter by about 1 mag dLaw et alj l2006l) . 
2M 1314+1320 has strong Ha emission, with an equiva- 
lent width of -54 A, corresponding to Lea/i'boi ~ -3.2 
dLepine et al.l2009h . consistent with the rotationally-saturated 
level observed in early- and mid-M dwarfs. 

' The VLA is operated by the National Radio Astronomy Observatory, a 
facility of the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 



2.2. VLA Observations 

We observed 2M 1314+1320 with the VLA on several oc- 
casions and at several frequencies (1.43, 4.86, 8.46 and 22.5 
GHz). Data were acquired as part of programs AB1245 and 
AB1312 (PI: Berger). The details of the observations are pro- 
vided in Table [T] All observations were obtained in the stan- 
dard continuum mode with 2 x 50 MHz contiguous bands. 
During the 10 hr simultaneous observation we used fourteen 
antennas at 4.86 GHz and twelve antennas at 8.46 GHz. Phase 
calibration was performed using J 1309+119, while the flux 
density scale was calibrated using 3C286 (J 1331+305). We 
reduced and analyzed the data using the Astronomical Im- 
age Processing System (AIPS). We inspected the visibility 
data for quality, and noisy points were removed. In each 
of the observations a source coincident with the position of 
2M 1314+1320 was detected. To search for source variability 
in the 10-hr simultaneous observation, we constructed light 
curves by plotting the real part of the complex visibilities at 
the position of 2M 1314+1320 as a function of time using the 
AIPS DFTPL routine. 

2M 1314+1320 is also identified in data from the Faint Im- 
ages of the Radio Sky at Twenty-Centimeters (FIRST) survey, 
where it was observed in December 1999 at 1.43 GHz with a 
flux density of 1 126 ± 123 /xJy. A comparison of the position 
in the FIRST data and in our subsequent VLA observations 
in 2009 leads to a proper motion of (5RA « -260 mas yr"' 
and (5DEC « -177 mas yr"', in excellent agreement with the 
known proper motion of 2M 1314+1 320. 

2.3. VLBI Observations 

We observed 2M 1314+1320 with the 10-antenna Very 
Long Baseline Array (VLBA) and the Green Bank Telescope 
(GBT) at 8.46 GHz on 2010 April 12 UT for approximately 
8 hr (program BM0327; PI: McLean). Phase-referenced ob- 
servations were carried out with a data rate of 512 Mbit s"' 
in dual polarization, using two-bit sampling. Eight base-band 
channels of 8 MHz bandwidth each were used, spread over at 
total base-width of 32 MHz. The correlator dump time was 
1 s. The complex gain calibrator was J 1309+1154, located 
at a distance of 1.8° from the target. There is a gap of about 
50 min in the resulting light curve of 2M 1314+1320 due to 
geodetic quasar observations intended to improve the astrom- 
etry. 

We detected 2M 1314+1320 with a flux density of /v = 
757 ±53 ^Jy at 8.46 GHz with 14±4% circular polarization. 
The source is unresolved with a beam size of 2 x 1 mas, cor- 
responding to an upper limit on the source diameter of about 
25 stellar radii (Figure [T]). This is about an order of magni- 
tude larger than the expected scale of the magnetic field. The 
flux density is lower by about 27% from the average of the 
VLA observations, most likely due to a loss of flux caused by 
phase calibration errors. We tested this hypothesis by elim- 
inating the longest baselines (i.e., reducing the angular reso- 
lution) and recovered the same flux density. In addition, we 
searched for potential emission from the secondary star in the 
system at a level corresponding to the missing flux (« 0.3 
mJy) but we did not detect any other source above a 4cr level 
of about 0.2 mJy within a 0.8" x 0.8" region centered on the 
detected source. 

2.4. Optical Spectroscopy 

To measure the projected rotation velocity, vsin/, we ob- 
tained a high-resolution spectrum with the Magellan Inamori 



Kyocera Echelle (MIKE) spectrograph mounted on the Mag- 
ellan/Clay 6.5-m telescope on 2009 May 19 UT. The sin- 
gle 1000 s spectrum was reduced using a custom reduction 
pipeline'* written in Python. Wavelength calibration was per- 
formed using ThAr arc lamps, and air-to-vacuum and helio- 
centric corrections were applied. We measured the rotation 
velocity by cross-correlating the sp ectrum with a templat e 
spectrum of the inactive star Gl 876 (iReiners & Basrill2008f) . 
The value of vsin/ was determined by comparing the width 
of the correlation function to cross-correlations between the 
template and artificially br oadened versions of th e template 
spectrum (see for example iReiners & Basrill2006h . We de- 
rive vsin/ from several spectral orders and find 45 ± 5 km s"' . 
There is no evidence for a spectroscopic binary. Ha emission 
was detected with an equivalent width of 9.9 A, correspond- 
ing to LjialLbox K, -A.Q, lower than previous measurements 
(lL6pine et al.ll2009h . 

We also observed 2M 1 3 14h- 1320 using the Low Dispersion 
Survey Spectrograph (LDSS3) on Magellan/Clay 6.5-m tele- 
scope on 2009 February 25 UT. A total of six 300 s integra- 
tions were obtained with the VPH-Blue grism using a 1" slit. 
The data were reduced using standard procedures in IRAF, 
and the wavelength calibration was performed using HeNeAr 
arc lamps. We detect non-variable Ha emission in all the in- 
dividual spectra with an average equivalent width of w 14.6 
A, corresponding to Lea/^boi ~ -3.8. 

2.5. Photometric Monitoring 

To study the photometric behavior of 2M1314H-1320 we 
observed the source on two consecutive nights (2010 April 
6 and 7 UT) using KeplerCam on the Fred Lawrence Whip- 
ple Observatory (FLWO) 1.2-m telescope for about 7 hr on 
each night. The observations alternated between the g- and 
/-band filters with individual exposures of 180 and 60 s, re- 
spectively. The data were reduced using standard procedures 
in IRAF. Relative photometry was performed usi ng the Dif- 
ference Imaging Photometry Pipeline photpipe (iRest et al.l 
1^05). 

Longer term photometric monitoring (2009 January 1 1 to 
2010 March 18 UT) was performed with a single telescope 
of the MEarth array at the Fred Lawrence Whipple Obser- 
vatory. MEarth consists of 8 identical 0.4-m robotic tele- 
scopes, and is dedicated to a search for super-Earth ex- 
oplanets in the habitable zones of nearby mid- to late-M 
dwarfs (iNutzma n & Charbonneau 2008). All observations 
were taken with a 715 nm long -pass filter, with the red end of 
the bandpass defined by the CCD quantum efficiency curve. 
Exposure times for data taken between 2009 January 1 1 and 
June 23 UT were 28 s with the telescope defocused to 6 pixels 
FWHM; data between 2009 December 5 and 2010 March 18 
UT were obtained with an exposure time of 1 1 s with the tele- 
scope operated in focus. Observations were scheduled at the 
default 20 min cadence until 2009 June 12 UT when the ca- 
dence was increased to 10 min. On 2010 February 14 UT, we 
observed for the entire night at the highest possible cadence 
(w 30 s, including overheads). 

Basic reductions and light curve generation were performed 
using an automated pipeline, based on the Monitor project 
pipeline (I rwin et an[2()07 ). A number of instrument-specific 
refinements have been made, and these will be described in 
full in a forthcoming publication (Berta et al. 201 1, in prep). 

'^ [http : //www ■ ociw . edu/Code/mike| 



MEarth data have two known systematic effects remaining af- 
ter the standard pipeline differential photometry corrections, 
which are of particular importance for the detection of con- 
tinuous photometric modul ations such as rotation, and are de- 
tailed in Ir win et al.l (120111) . These are the "meridian offset" 
and "common mode" effects. We follow the same method 
here to correct for these effects, fitting for them simultane- 
ously with the source modulations. 

3. PROPERTIES OF THE RADIO EMISSION 

The two initial 1-hr observations of 2M 1314H-1320 at 8.46 
GHz demonstrate a stable flux density level of about 1150 
/ijy and stable circular polarization of about 20%. Compared 
to the FIRST detection in 1999 (at 1.43 GHz), the flux density 
level appears to be stable for a period of about 10 yr. The 10- 
hr simultaneous observation obtained about 1 .5 yr later allows 
us to study the short-term temporal variability of the radio 
emission. The light curves are shown in Figure |2] The aver- 
age flux density is 1099 ± 18 ^liy at 4.86 GHz and 1032 ± 16 
/xJy at 8.46 GHz, again indicative of stable long-term radio 
emission, as well as a flat spectrum. The corresponding av- 
erage luminosities are L^(4.86) = (3.31 ±0.05) x 10"* and 
L^(8.46) = (3.52 ±0.06) x 10'"* erg cm'^ s"' Hz"'. Assum- 
ing that the emission comes from the primary object, the ratio 
relative to the bolometric luminosity is \og{vL^ / Lbo\) ~ -5.9. 
This value is about 5 times larger than the brightest radio- 
emitting M dwarfs detected to date, and roughly comp arable 
with the L3.5 dwarf 2M0036H- 1821 dBerger et alj|2005l) . 

The key feature of the 10 hr observation is the clear 
periodic variation in the flux density with P « 4 hr. 
Unlike the periodic radio emission detected in previous 
ultracool dwarfs (TVLM 5 13-46546 LSR1 8 35-H32, and 
2M0 746-H2000; iHalUnan et al] l2006l IIoOTI: iBergeretalJ 
120091) . the variation in flux density is smooth and sinusoidal 
rather than consisting of short-duration pulses. A Lomb- 
Scargle periodogram of the data is shown in Figure [3] re- 
vealing a single significant peak at both frequencies with 
P = 3.92± 0.06 hr (4.86 GHz) and 3.85 ±0.08 hr (8.46 GHz). 
This is consistent with a single period of 3.89 hr at both fre- 
quencies. Cross-correlation of the two light curves shows that 
there is no detectable phase shift between the two frequencies, 
with a 3(7 upper limit of « 3 min, or about 0.01 P (Figure [3]). 
We also fit the light curves with a sine-function using a least- 
squares algorithm, and find P = 4.0 ±0.1 hr (4.86 GHz) and 
3.9 ±0.1 hr (8.46 GHz). The ampHtude of the variation at 
4.86 GHz is 32 ±4% and at 8.46 GHz it is 22 ± 3%. We do 
not detect obvious periodicity during the 8 hr VLBA observa- 
tion, but note that the noise level in the light curve precludes 
the detection of variability at the level 20% as observed in the 
VLA observations. 

The periodic variability is also detected in circular polar- 
ization. At 4.86 GHz, the average polarization during the 
10 hr observation is close to zero but it varies at a level of 
24 ± 10% of the total flux density, alternating between right- 
and left-handed circular polarization over a best fit period 
of 3.81 ±0.09 hr from a Lomb-Scargle periodogram analy- 
sis (Figure H]). The periodicity of the circular polarization 
is correlated with the modulation of the total intensity, with 
the left-handed polarization coinciding with the peaks in to- 
tal intensity. At 8.46 GHz the average polarization is weaker, 
« 10%, with no significant periodicity, although the helicity 
of polarization may reverse during the peaks of the total in- 
tensity (Figure m. 

The radio spectral energy distribution (SED) between 1 .43 



and 22.5 GHz is shown in Figure |5] including all of our VLA 
observations and the FIRST survey data. Unfortunately, only 
short observations were made at 1.4 and 22.5 GHz, so we 
have no information about periodicity at these frequencies. If 
the flux density variations have a similar amplitude to that ob- 
served at 4.86 and 8.46 GHz, the measured flux may be biased 
low or high by about 20-30% depending on the phase of the 
light curve probed by these short observations. The broad- 
band spectrum peaks at about 5 GHz, with a flat spectrum at 
higher frequencies, /? sa -0.2 (/v oc ly^). For comparison, in 
Figure |6] we show the multi-frequenc y data for TVLM 513- 
46546 (M8.5) from 1.4 to 8.5 GHz JOsten et al.1 l2006h and 
DEN IS 1048-3956 (M8) from 4.9 to 24 GHz dRavi et alj 
l20Tlh . The broad-band SED of 2M1314H-1320 is signifi- 
cantly flatter than that of DENIS 1048-3956. 

The SED of circularly polarized flux is shown in Figure |7J 
including the full range from the 10 hr observation. It is un- 
clear if the measurement from the single 1 hr observation at 
1 .43 GHz represents the true mean polarization or if the ob- 
ject oscillates between right- to left-handed polarization at this 
frequency as it does at 4.86 GHz. Interestingly, the polariza- 
tion is consistently positive at 8.46 GHz. However, the 10 hr 
observation has the lowest average polarization, indicating a 
mixture of negative circular polarization related to the overall 
flux variability. The level of circular polarization in the VLBA 
detection (+14%) is consistent with the VLA observations. 

The radio emission properties can be used to extract infor- 
mation on the magnetic field strength in the context of the ex- 
pected emission mechanisms. The overall low degree of cir- 
cular polarization, and the wide frequency range over which 
the radio emission is detected may be indicative of gyrosyn- 
chrotron radiation. However, the nearly flat spectrum corre- 
sponds to an electron power law distribution with an energy 
index of p « 1.3, w hich is unus ually low compared to typical 
values of ;? ^ 3-4 (' Gudelll200 2). The magnetic field strength 
in the case of gyrosynchrotron emission can be estimated us- 
ing the relation f, w 2.85 x 10^B"■V;"■^ indicating B « 50 G 
for a peak frequency of j/^ « 5 GHz and a fraction of circular 
polarization of fc « 0.25. On the other hand, electron cy- 
clotron maser emission, which is produced primarily at the 
fundamental electron cyclotron frequency, i/f « 2.8 x lO^B 
Hz, require the magnetic field strength to vary from about 
0.5 to 8 kG within the emission region of 2M 1314+1320 to 
explain the broad-band SED. We note that the flat spectrum 
is indicative of an even wider distribution of magnetic field 
strengths. The maximal value is much larger than the aver- 
age surface magnetic field strengths measur ed on similar ob- 
jects using Zeeman measurements (1-4 kG: iReiners & Basril 
|2007|) . but it is possible that the Zeeman technique is not sen- 
sitive to such fields if their spatial scale is relatively small. 

4. PROPERTIES OF THE OPTICAL EMISSION 

The g- and /-band light curves from the FWLO 1 .2-m obser- 
vation are shown in Figure [8] The data exhibit the same peri- 
odicity seen in the radio, and the resulting phased light curves 
(with P « 3.78 hr) are shown in Figure|9] The variability am- 
plitude is « 80 mmag in g-band and « 50 mmag in /-band, 
and both bands are in phase. A phased light curve from the 
MEarth data is shown in Figure[TOl with a period of 3.785 hr 
determined from a Lomb-Scargle periodogram analysis (Fig- 
urefTTT). The variation amplitude in the 715 nm long-pass filter 
is w 15 mmag indicating an overall trend of decreasing vari- 
ability at longer wavelengths. 

The phased Ught curve is asymmetric, with a longer rise 



in brightness compared to the subsequent decline (the rise 
time from trough to peak is about 2.3 hr, while the decline 
back to minimum takes about 1.5 hr). The optical period is 
slightly shorter than the period determined from the VLA ob- 
servations, by 6 ± 3 min. This subtle difference may be po- 
tentially due to differential rotation if the radio and optical 
variability arise at different latitudes on 2M 1314+1320. The 
MEarth and VLA observations could not be phased due to the 
large gap in the observations. The overall sinusoidal variabil- 
ity of 2M 1314+1320 points to a large covering fraction for 
the spot(s), while the asymmetry is indicative of multiple or 
irregularly-shaped spots. 

Photometric rotational modulation is common in M dwarfs 
and is generally attributed to magnetic spots. However, 
the cool neutral atmospheres of ultracool are expected to 
also form clouds that may produce periodic variability 
(lBailer-Jonesll2002l) . Since 2M 1314+1320 is an M7 dwarf, 
its atmosphere is still fairly warm and dust clouds are less 
likely. The ratio of the variability amplitude in the g and 
/ bands is about 1.6, which is consistent with a cool spot 
(Rockenfelleret al. 2006). For comparison, the radio active 
M8.5 dwarf TVLM 5 13-46546 exhibits photometi-ic modula- 
tion with the same period as the radio variability, but the ^- 
and /-band variations are a nti-correlated, indicat ive of clouds 
rather than magnetic spots jLittlefair et alj|200 8). Here the g- 
and /-band variations are in phase, indicating that magnetic 
spots are the likely explanation. However, it remains unclear 
whether the optical and radio variability are produced by the 
same active regions or even related magnetic regions. 

5. TOPOLOGY OF THE RADIO EMISSION REGION 

With a rotation period of 3.89 hr and vsin/ = 45 ±5 km 
s"', the minimum radius of 2M 1314+1320 is inferred to be 
9 X 10^ cm (0.13 R^). Th e expected mass for an M7 dwarf 
iBaraffe & Chabrieii (il99g ) is about 0.09 M(r,, and the ex- 
pected radius is therefore about 0.1 Rq (Demorv et al. 2009). 
Thus, the inferred radius is larger than predicted, especially if 
/ ^ 90° , and this may be indicative of the trend for larger radii 
than model predictions (b y up to ^ 50%) for magnetically 
active ultracool dwarfs (Chabr ier et alj|2007l : iLopez-MoralesI 
2007; Jackson et al. 2009). Given that the minimum radius 
is already 40% larger than the model predictions, this sug- 
gests that that inclination is likely > 70°. However, since 
the binary companion is unresolved in /T-band observations, 
the magnitude of the primary is difficult to estimate. /T-band 
magnitude-mass relations indicate the mass may actually be 
higher (0.1-0.17 Mq) and the relatively low temperature may 
be a result of the inflated radius. This may resolve some of the 
discrepancy between the radius and Teff . 

The periodic radio emission, and especially the periodic he- 
licity reversals at 4.86 GHz, can be consistently explained 
with a simple geometric model: a large-scale dipolar mag- 
netic field mis-aligned relative to the rotation axis and with 
opposite polarity at each pole: see Figure [121 The inclination 
of the magnetic axis with respect to the observer leads to a 
larger projected surface area for one of the poles, and hence 
an increase in the total intensity when it rotates into view. At 
4.86 GHz we observe a helicity reversal when the other pole 
is in view. The fraction of circular polarization remains con- 
stant, but may actually be higher than observed if the inclina- 
tion is such that both spots always remain partially visible. If 
this is the case, the opposite polarities of the two spots would 
cause some of the circular polarization to cancel out. The 
lack of obvious helicity reversals at 8.46 GHz could be due 



to changes in the optical depth and the apparent size of the 
two magnetic poles as a function of frequency, particularly if 
higher frequencies are emitted near the magnetic foot-points. 
This is also supported by the smaller variability amplitude at 
8.46 GHz. 

The optical periodicity may be caused by cool spots spa- 
tially connected with the radio-emitting poles, or by unrelated 
spots, which are relatively common in M dwarfs. As noted 
above, spatial separation and differential rotation may explain 
th e small diffe rence in inferred optical and radio periods. 

iMorin et all (12010) propose two distinct magnetic topolo- 
gies of mid- to late-M dwarfs: large-scale, axisymmetric 
poloidal topology, and a complex, non-axisymmetric topol- 
ogy. Our geometric model indicates that 2M 1314H-1320 has 
a topology more similar to the large-scale poloidal fields 
rather than those with complex, small-scale fields . The reason 
for the se two populations remains unclear, but iMorin et all 
( 1201 Ol) propose that objects may transition towards larger- 
scale poloidal fiel ds with age . We note, however, that the 
objects studied in iMorin et aP (1201 0) have much slower ro- 
tation velocities (or correspondingly larger Rossby numbers) 
than2M1314H-1320. 

6. CONCLUSIONS 

2M 1314H-1320 provides a unique window into 
magnetically-induced radio emission in ultracool dwarfs, 
coinciding with the spec tral type at which the radio/X-ray 
correlation breaks down ( Bergerll2002l 120061) . It has the most 
luminous stable radio emission over the widest frequency 
range (1.4-22.5 GHz) to date, and furthermore exhibits clear 
sinusoidal periodicity, with the first example of a helicity 
reversal in the circularly polarized emission of an ultracool 
dwarf. Taken in conjunction, these properties point to 
emission from a simple dipolar field, mis-aligned relative to 
the stellar rotation axis. In the context of electron cyclotron 
maser emission, the magnetic field strength is required to 
be at least 8 kG. The inferred field structure is similar to 
the simple field topologies inferred fro m Zeeman Dopple r 
imaging of some fully convective dwarfs dMorin et alj|2010l) . 
while the field strength in the case of ECM emission is 
larger than the typical av erage fields measured from Zeeman 
broadening of FeH lines dReiners & BasriluOOTh . 

The stability of the radio emission for at least a decade in- 



dicates that the underlying dynamo is stable well beyon d the 
conve ctive turnover timescale for an M7 dwarf ('^ 70 d: ISaari 
1200 lb . Moreover, a comparison to convective dynamo models 
and simulations supports p redictioris of dominant large-scale , 
low-multipole fields (Dobleretal. 2006; Browning 2008), 
particularly for fast rotators (Browning 2008) . but is at odds 
with models that predict small-scale, high- multipole fields 
(iDurnev et al.lll993tlChabrier & Kukedl20()6b . 

Future observations of 2M1314H-I320 will address sev- 
eral key questions. First, if the emission is due to the ECM 
process, we expect a cut-off in the radio emission at a fre- 
quency that corresponds to the maximal field strength, j/max = 
2.8 X lO^Bmax Hz. Measurements of this possible cut-off fre- 
quency will therefore provide a crucial constraint on the mag- 
netic field since the Zeeman techniques are only sensitive to 
the average surface field or its large-scale component. The in- 
creased sensitivity and bandwidth of the EVLA, particularly 
at frequencies of ^ 20 - 40 GHz, is well-suited for this test. 
Second, simultaneous long-term monitoring of the radio and 
optical emission will test whether the variability has a similar 
period and phase. If identical, this will point to an origin in a 
related magnetic field structure, while a clear difference in pe- 
riod will be suggestive of differential rotation and an origin at 
different latitudes. Finally, long-term astrometric monitoring 
with the VLBA will be able to uncover a companion down to 
a few Jupiter masses (Forbrich & Berger 2009). Thus, further 
studies of 2M 1314H-1320 are bound to provide new insights 
into the properties of fully convective dynamos and the release 
of magnetic energy in the coolest stars and brown dwarfs. 
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TABLE 1 








VLA Observations of 2M 1314+1320 




UT Date 


Config. 


Duration 


Frequency 


Fu 


F^.v" 






(hr) 


(GHz) 


(MJy) 


(MJy) 


2007 Jul 20.02 


A 


1 


8.46 


1156 ±27 


215±26 


2007 Aug 10.00 


A 


1 


8.46 


1144 ±44 


227 ±43 


2009 Mar 25.12 


B 


10 


4.86 


1099 ±18 


38 ±17 


2009 Mar 25.12 


B 


10 


8.46 


1032 ±16 


117±16 


2009 Mar 27.32 


B 


1 


4.86 


1063 ±39 


<84 


2009 Mar 27.32 


B 


1 


8.46 


1052 ±42 


196 ±36 


2009 May 10.42 


B 


1 


8.46 


942 ± 30 


162 ±28 


2009 May 11.37 


B 


1 


1.43 


890 ± 58 


-218 ±6 


2009 May 11.41 


B 


1 


22.5 


763 ± 84 


<225 



Note. — " Circularly polarized flux density. 
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Fig. 1 . — Very Long Baseline Array (VLBA) 8.46 GHz image of 2M 1314+1320. The source is unresolved at the VLBA scale, with a beam size of 2 X 1 mas 
(bottom left). 
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Fig. 2. — Radio light curves for tlie 10 hr observation of 2M 1314+1320 at 8.46 GHz (top) and 4.86 GHz (bottom). The top panel at each frequency shows the 
total intensity (Stokes I) while the bottom panels display the circularly polarized flux (Stokes V). The 30-s averaged data are shown in grey, with 25-min averages 
shown as black circles. Sinusoidal periodicity is clearly visible at both frequencies. The best fit sinusoidal models are show in red, corresponding to a 22 ± 3% 
variation with a 3.9 ± 0. 1 hr period at 8.46 GHz, and 32 ± 4% variation with a 4.0 ±0.1 hr period at 4.86 GHz. The polarization hght curve at 8.46 GHz does 
not show significant periodic variation, while at 4.86 GHz it alternates between right- and left-handed circular polarization over a best fit period of 3.8 ± 0.4 hr 
with an amplitude of 24 ± 10%. 
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Fig. 3. — Lomb-Scargle periodograms for the total intensity data at 4.86 GHz (solid line) and 8.46 GHz (dashed line). A strong peak is seen at 3.92 ±0.06 hr 
and 3.85 ± 0.08 hr, respectively. Significance levels of 1 - lOcr are plotted as horizontal lines. Also shown is a cross-correlation of the data from both frequencies 
(inset), with a peak at 1 ± 1 min indicating no discernible lag between the two light curves. 
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Fig. 4. — Same as Figure[3]but for the circular polarization light curves. A strong peak is seen at 3.81 ± 0.09 hrs in the 4.86 GHz data, while no significant 
periodicity is evident at 8.46 GHz. Significance levels of 1 -8(t are plotted as horizontal lines. 
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Fig. 5. — Radio spectral energy distribution of 2M 1314+1320 from 1.4 to 22.5 GHz. Measurements from several short observations are shown in addition 
to the total flux density range detected in the 10 hour observations (vertical gray bars). The short observations at 1.4 and 22.5 GHz provide no indication of the 
variability and may not represent the mean intensity at these frequencies. The spectrum is surprisingly flat and appears to peak around 5 GHz with a spectral 
index of about —0.2 at higher frequencies. 
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Fig. 6. — Comparison of the radio spect ral energy distribution of 2M 1314+1320 with other ultracool dwarfs: TVLM 5 1 3-46546 nVI8.5: IOsten et al.l200^ and 
DENIS 1048-3956 (M8; (Ravi et aU2011l) . The SED of 2M 1314+1320 is significantly flatter than that of DENIS 1048-3956 
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Fig . 7 . — Same as Figure|5]but for the circularly polarized flux. It is unclear if the single 1-hr observation at 1 .4 GHz represents the true mean polarization or if 
the object varies from right- to left-handed circular polarization as it does at 4.86 GHz. The polarization is consistently positive at 8.46 GHz, but is lowest during 
the 10-hr observation. 
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Fig . 8 . — Optical g- and I'-band light curves from two nights of observations at the FLWO 1 .2-m. A period of about 4 hr is evident. 
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Fig. 9. — Phased light curves of the data in Figure [8]folded with a period of 3.78 hr. Data from the two consecutive nights are plotted as triangles and circles, 
respectively. The shape of the phased light curve appears more complex than a simple sinusoid. 
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Fig. 10. — A subset of high cadence photometric observations from the MEarth telescopes. A period of 3.78 hr is evident with an asymmetric light curve. 
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Fig. 11. — Lomb-Scargle periodogram of optical observations from the MEarth telescopes. A clear peak is seen at a period of P = 3.785 it 0.002 hr. Aliased 
or harmonic peaks are also visible due to the window function produced by the uneven sampling of the light curve (inset). 
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Fig. 12. — Simple geometric model of a mis-aligned dipolar magnetic field that can explain the various properties of the observed radio emission. The radio 
emission arises from two large magnetic hot spots with opposite polarities, located at opposite poles. The inclination with respect to the observer leads to a larger 
projected area from one of the poles (red) and hence a peak in the total intensity when it rotates into view. Helicity reversals are observed when the second pole 
(blue) is visible. 



